Electronic polarization at surfaces and thin films of organic molecular crystals: 

PTCDA 
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The electronic polarization energies, P = P+ + P- , of a PTCDA (perylenetetracarboxylic acid 
dianhydride) cation and anion in a crystalline thin film on a metallic substrate are computed and 
compared with measurements of the PTCDA transport gap on gold and silver. Both experiments 
and theory show that P is 500 meV larger in a PTCDA monolayer than in 50 A films. Electronic 
polarization in systems with surfaces and interfaces are obtained self-consistently in terms of charge 
redistribution within molecules. 



I. TRANSPORT GAP 

The electronic structure of organic molecular crystals is 
strikingly different from the conventional inorganic semi- 
conductors, such as Si, in that the electronic polarization 
of the dielectric medium by charge carriers constitutes 
a major effect, with energy scale greater than transfer 
integrals or temperature The transport gap E t 

for creating a separated electron-hole pair has a sub- 
stantial (1 — 2 eV) polarization energy contribution ||] 
and exceeds the optical gap by ~ 1 eV. Limited overlap 
rationalizes the modest mobilities of organic molecular 
solids. Devices such as light-emitting diodes, thin film 
transistors, or photovoltaic cells require charge transport 
and are consequently based on thin films, quite often de- 
posited on metallic substrates Organic electronics 
relies heavily on controlling films with monolayer preci- 
sion, on forming structures with several thin films, and on 
characterizing the interfaces. The positions of transport 
states and mechanisms for charge injection are among 
the outstanding issues for exploiting organic devices. We 
focus here on the electronic polarization of crystalline 
thin films near surfaces and interfaces. We find that 
electronic polarization, and hence E t , in a prototypical 
organic crystal is significantly different at a free surface, 
near a metal-organic interface, in thin organic layers, and 
in the bulk. 

Weak intermolecular forces characterize organic molec- 
ular crystals, whose electronic and vibrational spectra 
are readily related to gas-phase transitions [0J^]. Due 
to small transfer integrals, charge carriers are molecular 
ions embedded in the lattice of neutral molecules. The 
transport gap E t in the crystal is derived from the charge 
gap for electron transfer in the gas phase, 1(g) — A(g), 
which is the difference between the ionization potential 
and the electron affinity. But crystals have electrostatic 
interactions even in the limit of no overlap, and charge 
carriers are surrounded by self-consistent polarization 
clouds. In contrast to polaronic effects, electronic po- 
larization is instantaneous and directly affects the posi- 
tions of energy levels. Formation of polarization clouds is 



associated with stabilization energy P + for cations (the 
"holes") and P_ for anions (the "electrons"). The com- 
bination P = P + + P_ occurs in E t = 1(g) — A(g) — P. 
Since Coulomb interactions are long-ranged, polarization 
clouds extend over many lattice constants and P depends 
on the proximity to surfaces and interfaces. 
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Fig. 1 Schematic representation of charge-generation pro- 
cesses in crystalline molecular films, (a) UPS/IPES generates 
a cation/anion at the outer surface, while charge injection 
from the substrate involves the layer next to the metal, (b) 
Tunneling through a monolayer. Dashed ovals in (a) and (b) 
represent image charges in the metal. 

Figure 1(a) depicts schematically an ultraviolet pho- 
toemission (UPS) process where the ejected electron 
leaves behind a molecular cation in the outermost layer. 
In inverse photoelectron spectroscopy (IPES), the surface 
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is irradiated with low-energy electrons and the emitted 
photon is detected when an electron is captured to form 
a molecular anion. UPS data is increasingly available 
from sub-monolayer to ^100A films Q, while the IPES 
data is more limited. The combination of UPS with IPES 
yields E t directly, with P about 1 — 2 eV in representa- 
tive organic materials used in devices ||. As sketched 
in Fig. 1(b), tunneling electron spectroscopy gives E t as 
the interval between the differential conductance peaks 
when the potential of the tip matches either the electron 
or the hole transport levels of the film. 
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Fig. 2 (a) Composite UPS/IPES spectra as a function of 
PTCDA thickness on Ag. Energy scales are aligned by mea- 
suring the Fermi energy by UPS and IPES on Ag prior to 
PTCDA deposition, (b) dI/dV(V) STS spectra of filled and 
empty states recorded for a monolayer (bottom) and a 2-3 
molecular layer (top) film of PTCDA deposited on Au. The 
corresponding STM images of the films are shown. The curves 
were recorded at the same tunneling setpoints as the corre- 
sponding area scans. 

Figure 2(a) shows UPS and IPES spectra of PTCDA 
(perylenetetracarboxylic acid dianhydride) on silver. 
PTCDA is an excellent former of crystalline films whose 
structures are close, though not identical, to having a 
(102) plane of the bulk crystal in contact with the sub- 



strate j|. The measured transport gap on thick (> 50A) 
films is Ef = 3.8 eV on Ag, in excellent agreement with 
PTCDA on Au |. We use PTCDA on Ag rather than 
Au because the Au(5d) levels interfere with UPS of a 
monolayer. Ef includes a 0.2 eV intramolecular vibra- 
tional contribution [|| that reduces the UPS/IPES gap. 
This correction is the same for all PTCDA films. Care- 
ful analysis of peak positions indicates ~200 meV shifts 
of both the cation level (UPS) and anion level (IPES) 
between mono- and multilayer films. 

Figure 2(b) shows the results of scanning tunneling 
spectroscopy (STS) on a monolayer (lower spectrum) and 
a 2-3 layer film (upper spectrum) of PTCDA on Au(lll). 
These dI/dV(V) spectra represent the density of filled 
and empty states involved in tunneling out of and into 
the layer, respectively. Each spectrum is the average of 
25 spectra recorded at various points on highly ordered 
molecular layers. High-resolution scanning tunneling mi- 
croscopy (STM) images of these layers, taken concomi- 
tantly with STS, show the characteristic difference Q] 
between monolayers, in which molecules are in contact 
with, and parallel to, the Au surface and appear sym- 
metric, and second and subsequent layers, in which the 
tilt angle of molecules introduces an asymmetry in their 
STM image. 

The tunneling spectrum of the monolayer shows peaks 
leading to Ef 11 " — 3.3 eV for the energy difference be- 
tween adding an electron or hole. Remarkably, this en- 
ergy difference increases by about 0.25 eV on the 2-3 
layer spectrum, in excellent agreement with the interme- 
diate UPS/IPES spectra of a 16A film (Fig. 2a). Each 
peak shifts away from the Fermi level by a roughly equal 
amount with increasing coverage. We note that Ef and 
Ef — E^ show no dependence on dipoles at the metal- 
organic interface, which have opposite signs for PTCDA 
on Ag and Au g. The dipole at the PTCDA/ Ag inter- 
face corresponds to electron transfer from the metal to in- 
terface molecules, which gives rise to the filled gap states 
at -0.6 and -1.8 eV on the 4A UPS spectrum. The 
broad feature around -0.6 eV on the 4A PTCDA/Au 
STS spectrum corresponds to the Au surface state, which 
is not eliminated by the deposition of organic molecules 

Reduced polarization energy at surfaces has long been 
appreciated on general grounds jl],|| . As anticipated and 
found for thin films of anthracene on Au, P+ is about 
200 meV smaller in the surface layer [|lo]jll]]. With simi- 
lar reductions expected for P_, the transport gap at the 
surface is increased some 400 meV from its bulk value, 
E t . The gap increases near surfaces because vacuum is 
not polarizable. Conversely, E t decreases near organic- 
metal interfaces due to the high polarizability of metals. 
These opposite contributions partially cancel in very thin 
films on metal substrates. 
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II. CHARGE REDISTRIBUTION 

We note that E t , which has been the focus of theoreti- 
cal study Jl2|,|l| and governs bulk transport, is beyond the 
reach of the surface experimental techniques that access 
the outermost layers of the material and thus reflect Ef . 
Similarly, the transport gap Ef 1 of the layer next to the 
metal is relevant for charge injection (Fig. 1(a)). A gap 
reduction of several hundred meV is significant since it 
directly alters interface barriers. 

No theoretical treatment exists for the electronic po- 
larization near surfaces and interfaces, which require an 
accuracy of ~100 meV. Methods to estimate P± in the 
bulk have been developed, primarily for the acenes, based 
on the microelectrostatics of polarizable points that rep- 
resent organic molecules 
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Dielectric response of 
a neutral surface has also been studied [0. We have 
recently developed an approach based on the analysis 
of charge redistribution in organic molecules |ll|, and 
demonstrated accurate calculations of P+ , P_ , and ener- 
gies of ion pairs in bulk PTCDA and anthracene crystals 
p6| . We apply here the same approach to calculate po- 
larization in thin organic films. 

In the zero-overlap limit, molecules comprising the 
crystal are quantum-mechanical objects interacting by 
classical forces. A self-consistent problem can be formu- 
lated that treats molecules rigorously in the external 
fields of all other molecules. We describe charge redistri- 
bution in organic molecules in terms of the atom-atom 
polarizability tensor Ely that relates a change in the par- 
tial charge at an atom i due to the electrostatic potential 



(rj) at atom j: 



plane at z = 0, a distance h from the innermost molecu- 
lar layer. Image charges and dipoles at z < ensure that 
the potential cj)(z = 0) = 0. Any r/>(0) = C leads to the 
same result for P = P + + _P_, since one charge is sta- 
bilized and the other is destabilized. The metal-organic 
separation h is a model parameter that is related to Van 
der Waals radii. 
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Fig. 3 Electronic polarization P — P+ + P_ at the outer 
layer of an JV-layer film (Fig. 1) at separation h from the 
metal. The P values at the free surface (N — > oo), a mono- 
layer (N = 1), free-standing films (h — * oo), and the bulk are 
indicated. A pill-box geometry with R ~ lOOA is used. The 
inset shows the transport gap, E\ n \ across a 10-layer film 
with n = 1 next to the metal and n = 10 at the outer surface. 



pf^ are the atomic charges in an isolated molecule. The 
tensor Hy is a natural extension of the similar quantity 
in 7r-electron theory jl7j. We compute Ily using INDO/S 
p8[ , which is a semiempirical Hamiltonian designed for 
spectroscopic molecular properties. Ily describes the ma- 
jor, "charge-induced" part of molecular polarizability a , 
which is augmented to reflect the actual polarizability a 
by introducing induced atomic dipoles \ii and distribut- 
ing the difference a = a — a c over 38 atoms of PTCDA in 
the spirit of submolecular methods Self-consistent 
equations for pi and fii are then solved for increasing 
cluster sizes of mesoscopic dimensions (^100A), and the 
macroscopic limits are found. For neutral lattices the ap- 
proach has yielded accurate anisotropic dielectric tensors 
of two representative organic crystals Jl5|] . 

We use identical molecular inputs here to model 
PTCDA films as infinite slabs terminated by (102) planes 
next to a vacuum and a metal, respectively, as sketched 
in Fig. 1. The metal is taken as a constant-potential 



We use a pill-box geometry, defined by the radius R 
(sketched in Fig. 3) with the ion placed in one of the 
layers near the center of the box. We initially find self- 
consistent atomic charges and induced dipoles ~p\ and Jif , 
k = 1, N in a neutral film of N layers, where transla- 
tional symmetry gives rapid convergence. We then con- 
sider pill-boxes of increasing diameter 2R to find self- 
consistent pf — ~p\ and [if — Ji!? for every molecule a in 
the pill-box. The largest systems (2i? = 135 A , N = 10) 
contain ^2400 PTCDA molecules and their images, and 
larger values of R can be used. P + is the energy differ- 
ence between a neutral film and one containing a cation. 
P_ is the corresponding energy for an anion in the same 
position. We note that the energy of an infinite film is 
extensive, but the difference is finite and can be evalu- 
ated ]l6[ ]. All data reported below is in the limit R — > oo. 
P thus depends on which layer contains the ion. 
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III. RESULTS 

For the ion in the outermost molecular layer in the 
limit N — > oo we obtain the free (102) surface polariza- 
tion P s = 1.41 eV. This limit does not depend on h or the 
metal. P s is 0.41 eV less than the bulk value P = 1.82 
eV Jl6[ . The difference is consistent with experimental 
estimates [M and corresponds to the surface correction 
c = 1 — P°/P = 0.23, where the value inferred || from 
UPS and IPES spectra was c - 0.25. 

Figure 3 gives results for PTCDA layers of finite thick- 
ness. Unlike bulk or free-surface calculations, which are 
essentially parameter-free, the finite layer data depends 
on h, which is the only way the metal enters our ideal- 
ized model. The reasonable value h — a places the metal 
plane one lattice spacing a = 3.214 A from the innermost 
layer, which is also consistent with Van der Waals radii, 
and 10 % variations of h are shown in Fig. 3. We also 
computed P at h — 1.244a, 2a and 4a, and extrapolated 
as 1/h. The limit h — > oo gives the polarization at the 
surface of a free-standing film of N layers. 

All curves in Fig. 3 converge to the free-surface value 
P s = 1.41 eV. Using the curve h — a we find the sin- 
gle monolayer value p ML = 1.93 eV, which corresponds 
to the tunneling spectroscopy setup (Fig. 1(b)). We see 
that for a monolayer on the metal surface, the polariza- 
tion energy is indeed close to the bulk value in line with 
the expected cancellation discussed above. The differ- 
ence P ML — P s = 0.52 eV agrees with the experimental 
Ef - £ t ML = 0.45—0.50 eV for PTCDA on Ag or Au. 
In fact, the agreement is slightly better since UPS/IPES 
data is for the the films of finite thickness N ~ 20. The 
free-standing monolayer has P — 0.64 eV, which is about 
Et/3. Such a big polarization is consistent with large in- 
plane polarizability of PTCDA molecules. These results 
for P are summarized for comparison in Table I. 

Table I. P+ + P_ at various positions in the film 

Bulk 1.82 eV 

Free surface 1.41 eV 

Monolayer, free-standing 0.64 eV 

Monolayer on metal (h = a) 1.93 eV 

Surface of a bilayer on metal (h = a) 1.73 eV 

Layer next to metal, thick film (h = a) 2.21 eV 



Analysis of UPS data for films of the electron-transport 
molecule Alq3 [tris(8-hydroxy-quinoline)-aluminum] on 
silver Jisj ] yields strikingly similar conclusions: the trans- 
port gap for a monolayer on metal substrate, P t ML , is 
equal to the inferred bulk value, and is about 400 meV 
narrower than the gap in the outermost surface layer of 
a thick film, Ef. Also, the inferred gap in the innermost 
layer of a thick film on metal surface, Ef 1 , is about 400 




meV less M than E. 



ML 



1 23456789 10 

Layer number 

Fig. 4 Variation of the transport gap from the bulk value, 

In) 

E t — Et, across a 10-layer film with n — 1 next to the metal 
and n — 10 at the outer surface. 

Figure 4 shows the variation of the transport gap 
E t across a 10-layer PTCDA film on a metal substrate 
(h = a). We note that surface effects extend several 
layers into the sample. The long-range nature of sur- 
face polarization has been largely neglected. An influen- 
tial early UPS study of 20A vapor-deposited anthracene 
films on Au ascribed the 200 meV shift of P+ to the sin- 
gle outermost layer j[o|. The additivity of polarization 
contributions suggested Jl9| for Alq 3 is tacitly based on 
short-range interactions. Greater polarizability next to 
the metal is consistent with strong charge confinement to 
interfaces, as inferred recently for pentacene field-effect 
transistors with remarkable electronic characteristics [g| . 
In general, the 400 meV increase of P+ +P_ at the metal- 
organic interface is not shared equally by the electron and 
hole. The stabilization of either carrier by roughly 200 
meV is important for matching transport levels in injec- 
tion. 

Weak overlap in molecular solids constitutes, in fact, 
a significant simplification over covalent bonding in in- 
organic semiconductors. Even though charge transport 
is critical for electronic applications, the picture of local- 
ized carriers is the proper zeroth-order approximation, 
to which the overlap (i.e. kinetic energy) should be con- 
sidered as a perturbation. This docs not by itself rule 
out the band-like description. Rather, charged quasipar- 
ticles are to be understood as surrounded by polarization 
clouds, likely affecting the band widths. 

Our results for P + and P_ are exclusively electronic. 
Lattice relaxation around charges are considered to be 
small corrections (~10%) on general grounds pL0]. The 
idealized model of image charges does not depend on the 
metal's Fermi energy or on surface dipoles, whose shifts 
cancel in P = P + + P_ . It also ignores surface states 
or surface reactions that are known to occur at specific 
organic/metal interfaces p] . Since the self-consistent cal- 
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culation requires the film's structure, it is not directly 
applicable to amorphous or structurally uncharacterized 
films. 

Support of this work by the NSF (DMR-0097133) is 
gratefully acknowledged. 
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